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We present very simple non-integral relations between deep inelastic structure func- 
tions Fjj, F2 and the gluon distribution at small x based on perturbative QCD 
which are useful for the phenomenological analysis of data at low x. As an appli- 
cation we extract the deep inelastic scattering cross-sections ratio R = ol/&t hi 
the range 10 -4 < x < 10~ 2 from F 2 HERA data. 



1 Introduction 

For experimental studies of high energy hadron-hadron and lepton-hadron pro- 
cesses it is necessary to know in detail the values of the parton (quark and 
gluon) distributions (PD) of nucleons, especially at small values of x. Of great 
relevance is the determination of the gluon density at low x, where gluons are 
expected to be dominant. 

The basic information on the gluon structure of nucleons is extracted from 
the measurement of the deep inelastic structure function F2 in lepton-hadron 
scattering (DIS). One of the usual procedures compares experimental data 
with the theoretical prediction for F2 obtained from the solution of a system of 
complicated coupled integro-differential quark and gluon evolution equations. 
It is also possible to extract the gluon distribution more directly from r Fi.scaling 
violations using a very simple relation with the Q 2 derivative of F2 lTEI. 

By other part, Fl or the ratio R — Fl/(F2 — Fl,), is also an interesting 
quantity, because it is a very sensitive QCD characteristic. For examplle, future 
Fl measurements^], will be used as a signal of the gluon structure at low xU. 

In perturbative QCD, there is the possibility to connect Fl with F2 due to 
the fact that at small x the DIS structure functions depend really on only two 
independent functions, the gluon and the singlet quark distribution (the non- 



a At present there are only preliminary measurements by HlE 



singlet quark density is negligible at small x), which in turn can be expressed 
in terms of Fi and its derivative dFijdlnQ 2 . 

In this article we present very simple linear relations || between the gluon 
density and F L (x, Q 2 ) with F 2 (x, Q 2 ) and dF 2 (x, Q 2 )/dlnQ 2 at small x. Using 
these formulas we exploit the possibility of extracting information about the 
gluon distribution and Fl at small x, directly from the measurement of the 
Fi scaling violations. This method complement the standard analysis where 
quarks and gluons, determined from complex fits to data, are integrated for 
the calculation of Fl ■ With our formulas it is possible to take into account the 
experimental uncertainty in the theoretical calculation more directly. 

The standard intial form for the singlet quar-k-sfx, Q 2 ) and gluon g(x, Q 2 ) 
distributions |] at some Q\ are parameterized bylliJ: 

p{x,Ql) = A v x- s >{l-xy*(l + e p y/x + *i P x) {p = S,g) (1) 



Until the recent tuoA the value of S p was a matter of discussionlHrEj, but the 
new HERA dataliil'E-3 start to overcome this controversy. From the theoretical 
side, the type of esjokdjion of the PD in Eq. [I] depends on the value and form of 
S p . For example ESO'El, a (^-independent 5s = S g obey the DGLAP equation 
when x~ Sp 3> 1. However, if S P (Q 2 ) — in some point Q 2 . > lGeV 2 , the 
behaviour p(x, Q 2 ) ~ Const is not compatible with DGLAP evolution and a 
more singular behaviour is generated HEJO. These cases have been recently 
evolved to a common picture Ej where partons are really a combination of 
two solutions (at x~ Sp 1 and at S p ~ 0) linked at some Q 2 point. For a 
Regge-like form of structure functions, one obtains0p(x, Q 2 ) ~ a; _ *i>W a ) with 
next-to-leading order (NLO) 8s{Q 2 ) ^ $ g (Q 2 ) intercept trajectories. Without 
any restrictions, it generates the double-logarithmical behaviour, 



p(x,Q 2 )~exp[2JS p (Q 2 )ln~) (2) 



where at NLO and for / = 4 active quarks one hasEl: 
being t = ln(a(Q 2 )/a{Q 2 )) and r = a{Q 2 )-a(Q 2 )$ 



'Some of them have been already published^ 
c We use PD multiplied by x and neglect the nonsinglet quark distribution at small x. 
^Hereafter we use a(Q 2 ) = a 8 (Q 2 )/4flr. Sometimes it is denoted as a to shorten long 
formulae. 



2 



2 The gluon and Fl as a function of F 2 and dF 2 /dlnQ 2 



By lack of space we will only presen 
the details will be given elsewhere E 



the final formulas of the calculation while 
Assuming the Regge-like behaviour of 



Eq. |l| for x~ Sp ^ 1, one can replace at small x the convolution integrals by 
ordinary products in the Q 2 evolution equations and in the formulas relating 
F 2 and Fl to PD. From the F 2 equation one can extract the singlet quark 
combination S(x,Q 2 ) as a function of F 2 (x,Q 2 ) and g(x,Q 2 ) and substitute 
it into the equations for Fl and for dF 2 /dlnQ 2 . 

The case of the non-Regge type behaviour has to be treated independently, 
but it is possible to combine both cases (Regge and non-Regge types) in a single 
formula valid for any value of S p : 



g(x,Q 2 ) 



2/ 



Isg 



~(1),1+<5 S 
■Isg 9 a 



0{a 2 



„l-<5x 



F L (x,Q 2 ) =aB s L M F 2 {x,Q 2 )-2 



a 



F 2 (x,Q 2 ) 



dF 2 {x,Q 2 ) 
dlnQ 2 



^ 1+Ss F 2 (x,Q 2 ) 



aR 



(0),l+<5, , ~(1),1+S g 

+ 7s 9 a 



^Sg 

+ 0{a 2 .x l - s ?) 



dF 2 {x,Q 2 ) 
dlnQ 2 



(3) 



(4) 



where e = e i is the sum of squares of / quark charges. The variables B^ v 
(k = 2,L) and {p,l = g,S) are respectively the one loop parts of the 

Wilson coefficients and anomalous dimensions of the operators. The second 
order quantities arc related to the two loop Wilson coefficients and anomalous 
dimensions by: 



H L - H L 

_(1) (1)4+5, 

TSg = TSg 



B, 



9,1+5, D S,l+<5 



B 



S,5 S (0),l+<5, 



'jB 



7 S 



B. 



9A+S S 



(2/3, 



o + 7. 



(0),l+«5, 



All these coefficients have to be continued analytically from the well known 
integer values n to the non-integer ones rj. 



The values of 7g g and R 9 L coincide with j Sg and R L , respectively, with 



the replacement, 



l dy g(y,Q 2 ) 

x v g{x,Q 2 ) 



(5) 
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Figure 1: The ratio R extracted from HI data. The lines are the predictions from different 
models (see text). The band represent the uncertainty from a DGLAP analysis of HERA 

data 



In the cases x Sg Const and 8 g — > 0, the r.h.s. of Eq. [| leads to l/<5 9 and 
1/Sg, respectively, where 



lln(l/x) 



S P (Q 2 ) A6 P {Q 2 ) 



O 



8Ji6 p (Q 2 )Hl/x) 



(_ 1 ) 

V4(5p hx{l/x)J 



(6) 



The comparison of our results with others is performed elsewhere E 



3 The extraction of the ratio R at low x 

As. an application we have extracted the ratio R(x,Q 2 ) frorp-LU 1994 
datali-3, determining the slopes dF2/dlnQ 2 from straight line-fits BE-a The 
systematic errors have been taken from an early analysis of HlEil. The running 
coupling constant a s (Q 2 ) has been calculated at two loops using A" (MS) = 
225MeV. 

Figure 1 shows the extracted ratio R at Q 2 — 20 GeV 2 for 6s = 0.3 and two 
different exponent for the gluon, 6 g = 6s and 8 g = 6s + 0.05. These value 
very close to those obtained by various groups from QCD fits to HI datalii 
Fig. 1 also shows some experimental data pointsES at high x. For comparison 
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we have also plotted various predictions for R using QCD formulae at 0(al) f| 
and parton densities extracted from fits to HERA data. The large difference 
between the result from MRS(G) and the latest set MRS(Rl) E3 shows the 
sensitivity of R to the update of these parton densities to new HERA data. 
One can notice that there is good agreement between MRS(Rl) prediction and 
our points when S g = Ss + 0.05 = 0.35, even though it is not surprising because 
these values were taken from the MRS analysis lj. However our result gives 
extra information about the uncertainty. A more precise future measurements 
at low x should lie within the error bars of the results presented in Fig. 1. 

By other part recent theoretical predictions on R b^ed on conventional 
NLO DGLAP evolution analysis of 1 
picture of BFKL dynamics (NPRW) I 



IRA data (LBY) Ea and on the dipole 
both finding values 5s = 5 g m 0.3, lie 



between the two of our above cases. 
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